Coulomb interaction and first order superconductor-insulator transition 
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The superconductor-insulator transition (SIT) in regular arrays of Josephson junctions is studied 
at low temperatures. Near the transition a Ginzburg-Landau type action containing the imaginary 
time is derived. The new feature of this action is that it contains a gauge field "I> describing the 
Coulomb interaction and changing the standard critical behavior. The solution of renormalization 
group (RG) equations derived at zero temperature T = in the space dimensionality d = 3 shows 
that the SIT is always of the first order. At finite temperatures, a tricritical point separates the 
lines of the first and second order phase transitions. The same conclusion holds for d = 2 if the 
mutual capacitance is larger than the distance between junctions. 

PACS numbers: 74.40.Kb, 74.81.Fa, 74.25.Dw, 64.60.Kw 
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Introduction -After several decades of intensive stud- 
ies Josephson junctions arrays (JJA) still remain an im- 
portant and inspiring problem. These, at first glance, 
simple systems exhibit, depending on parameters, super- 
conducting, insulating or metallic properties (see for a 
review [^). As a model, JJA is relevant for granular su- 
perconductors and disordered superconducting films [2]. 

Coulomb interaction (CI) is crucial for properties of the 
JJA. It suppresses the density fluctuations, i.e. causes 
fluctuations of phase (p of the superconducting order pa- 
rameter, which can destroy the superconductivity even 
at zero temperature, T — Q. As the charging energy of 
one grain increases, the JJA undergoes a superconductor- 
insulator transition (SIT) - an example of a quantum 
phase transition. As we show in this Letter the long- 
range nature of the CI qualitatively affects the SIT. 

JJA can be described by the effective Hamiltonian f3^ 



H = - N^ B\j._j.i\hrhr' — 2, Jr,r' COS {ifi^. 



Vr' 



(la) 
where Jr,r' is the Josephson energy of the junction be- 
tween neighboring grains r and r', and h^ is the particle 
number operator fij. = —idjdkp^ in a grain. We assume 
that superconducting pairing within each grain is strong 
and consider the limit of the inflnite single-electron gap. 
Accordingly, the excitations of isolated grains have the 
charge -|-2e/(— 2e) (bosons/antibosons). Matrix Br-r> 
describes the interaction between two bosons on the 
metallic grains r, r'. It is well approximated as 



B: 
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Here t is the period of JJA, and C > ^ is the mutual 
capacitance of the neighboring grains. The energy to 
charge one grain can be estimated as B^ ~ Br^e- 



Equation (la) looks like a Hamiltonian of a quantum 
XY -model whose critical behavior is described by an 
N = 2 -component (j) field theory. Such a field theory 



was discussed, e.g., in Refs. |4H6] for both finite T and 
T — 0. For T > 0, the critical behavior near the phase 
transition is described by a d-dimensional 2-component 
(f) theory: the quantum transition at T = involves 
imaginary time r as an additional dimension, i.e. the 
same ((> theory should be considered in (i+ 1 dimensions. 
In both cases the SIT is of the second order with the 
critical behavior of XY model in d or rf + 1 dimensions. 



This is correct if the matrices Br-r' in Eq- ( lb ) are ei 



j(0) 



ther diagonal, B^ /, — BoSry or sufficiently short ranged. 



The effect of the long range part of Eq. ( lb ) on the crit 



ical behavior has not been investigated so far. 

In this Letter we derive a field theory that properly 
describes the SIT at low temperatures in the disorder-free 
model (fTl). The Coulomb interaction (1/ \r — r'\ decay of 
Br,r' at large distances) results in an additional gauge 
field $ in the Ginzburg-Landau (GL) expansion near SIT 
and, at T — > 0, causes additional logarithmic divergences 
in the upper critical dimensionality d — 3 of JJA. We 
derive and solve renormalization group (RG) equations. 
Solutions demonstrate the first order SIT at sufficiently 
low temperatures. To describe SIT in 2D JJA we use 
e-expansion at small e = 3 — d. We find first order SIT 
as long C S> £, and it may become continuous otherwise. 

Problems with mean field description- To derive the 
field theory of fiuctuations near the phase transitions one 
should write a GL expansion in superconducting order 
parameter A = J (exp (i^?^)) , where (...) is the quan- 



tum mechanical average with the Hamiltonian (la), and 
J = ^j, Jr^r' ■ The mean field approximation is obtained 
by minimizing this expansion. At first glance, the free 
energy functional F [A] can be derived straightforwardly 
and should have a form of the standard GL expansion. 
Indeed, F [A] for time-independent A^ has the form 



F[A] 



r,r' 



\aryl^rK'+(irr'\^rV\^r 



(2a) 



where the phase transition is controlled by 

ar,r' = [j'^]^y-25r,r'/Eo, (2b) 

and i?o == Bq/2 is the energy of adding or subtracting one 
Cooper pair (boson) to a grain. Apparently, according 



Eqs. (2a|-(|2b|, there exists the critical coupling Jc 



Jc - So/2, 



(2c) 



so that at J > Jc the system is a superconductor, while 
the insulating state corresponds to J < Jc- However, this 
conclusion relies on a quartic term being local and pos- 
itive. Explicit calculation of the function (3^. ^, starting 



from Eq. (la) gives [7 
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(2d) 



where E^^, — Bq ± Bry are the energies of two bosons 
(boson-antiboson) pairs located on the grains r and r' . 
At large distances Br.r' ^ Bq, and Eq. (2d I yields 



/?!. 



-r'|>£ 



-bI^/eI, 



(2e) 



As at large distances B^-y ex |r — r'\ , the sum over r' 
in Eq. ( |2a[ ) is negative and diverges linearly (logarithmi- 
cally) for coordinate-independent |A| in 3(2)-dimensional 
JJAs. Such divergences at large distances signal that 
the additional soft modes should be included to make 
the theory local, and the conventional naive mean-field 
is not conclusive. This problem arises at low tempera- 
tures, T <C -Eo, only, while for T 3> i?o, the Coulomb 
interaction can be neglected and the SIT temperature in 
the mean-field approximation [3j turns out to be 



Tc = J. 



(3) 



In the vicinity of T^ the conventional GL free energy with 
a time-independent A is valid. 

Effective field theory- To account for the long range 
interaction we modify the model slightly: we separate 
Br,r' Eq- ( lb ) into the local and long-range parts and 



smoothen the latter: 



Br — 2Eo5r.O + B\j.\+jl, 



(4) 



where 7^1 which controls the short distance cut-off 
that will drop out of final results. Although we assume 
Eq 3> £drB\j,\\r~e, the number of long range terms is 
infinite and their effect accumulated from large distances 
has to be included together with quantum fluctuations of 
A. We separate the Hamiltonian into the bare one, Hq, 
and perturbations Hjb 



H — Hq + Hj + Hb, Hq — 2_^ EqU^, 



(5) 



hb=j: 



- , Hj = -'^ Jr,r' COS {pr- 'Pr')^ 



and write the partition function 

Z^Tr e-"/^ZQ ^T^e'/o^'" ''^[^■^(^)+^-M] ) . (6) 

Here T-r stands for imaginary time ordering, (. . • )o = 
Zo'^Trie-^°/T...\, Zq = Trle-^'^lrX, and 

Bb,.] (t) = e^«^HB,.je'"°''. Terms Hb (r) and Hj (t) 
are decoupled by Hubbard-Stratonovich fields A,, (r) 
(complex) and $r (t) (real): 

— = fvArVAlV'i're-'^^ (Tre-J'o'^ "^'^A , ^'^^ 

Zq J \ /o 

Hi=Y.{' ^f'r'^r] (r) + [Are-^^^] (r) + [A^e^*-] (r) }, 

r 

^1 = E /' d-" {-^^T^' ^^rK'] (r) + Ib;), [$.$,,] (t)| 

where J^r' KyB\r'\+ii = '^'■■O- 

In the limit of low temperatures T <^ Eq the average 
{■■■)q can be calculated over the ground state when all 
grains are neutral. We neglect exponentially small con- 
tributions like exp {—Ei/T) , where Ei are the eigenen- 
ergies of Hamiltonian Hq corresponding to the charged 
isolated grains. At the same time, we keep finite 1/T 
when integrating over r in Eqs. (l7| to obtain algebraic 
in T contributions. This calculation is carried out near 
the SIT by cumulant expansion assuming that the fields 
A (r, r) , ^(t, r) are slow in both time t and coordinate r. 
As we compute averages with the bare Hamiltonian Hq, 
all generated terms remain local. Introducing continuous 
coordinate description we obtain 



Z:^ I e- ^'>'^ dT{L^+L^)jy^^^^ ^)p^* (^^ T)V^{r, t). 



The order parameter is controlled by Lagrangian 
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(8a) 



(8b) 



where e — 3 — d. The coefficients in this action are ex- 
pressed in terms of the initial constants of the Hamilto- 
nian ^ at \dJ\ <^ Jc = EqI2] SJ = J - Jc, as 



B2/2 
2 ^0* 

c = -^ — ; a 
2d ' 



2SJ 



.f^=^.b=m^ (8c) 



r^r' 



Energy E is the running high-frequency cut-off in the 
theory, it starts at E = Eq. We included this cut-off 
explicitly in the action to keep the interaction constants 
dimensionless, avoid rescaling of (r, r) during RG, and 
explicitly illuminate the dimensionality of the interaction 
terms. As / can be always removed by the rescaling of A, 
the physical results may depend only on constants a, b, c. 



The fluctuating voltage $/ (2e) is of great importance 
for the critical behavior. This field is controlled by the 
Gaussian Lagrangian 



1 



drdr'^ (r, r) 



B- 



^{T,r'): 



(8d) 



After the coordinate Fourier transform of the fields 



^(k.r) = $*(—/(;, r), we obtain from Eq. (lb) 



(27r) 
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Harmonics ^{k) with |fc| > l/{ji) are suppressed. The 
coupling constants in Eq. ( 8e ) are defined as 



gf 



V6e 



d=3 n d=2 

92 = 0; 92 



gf 



n^EoC 



ttIEc 



(8f) 



For this approach to be applicable, the fluctuations of <& 
should be small i.e. 5i + 52 ^ 1- 

The field theory ^ defined for two slow fields A,<i> 
helps one to avoid the negative and non-local quartic 



terms Eq. (2e), arising in a single field formulation (2a I 
Indeed, Eqs. (I8|) are invariant under the transformation 



$^(r)^$^(r) + 0(T); A ^ Ae^'-^'f^ '^^^'^(^i) (9) 

for J„ (j)[T)dT ~ (mod 2tt), i.e. the effect of the fiuc- 
tuations of ^^e''^'^ at kc ^ u vanishes. Fixing by hand 
A while allowing for all the fiuctuations of $ violates the 
gauge invariance ([9]), and overestimates the contributions 
of small k and leads to an incorrect non-local theory. 

Renormalization group analysis- To the best of our 
knowledge, the model (Is]) has never been discussed. We 
analyze it using the RG approach in d = 3 and d = 3 — e. 



The cut-off dependence, Eq. in Eqs. (8b I, (8e| sug- 
gests that the theory is logarithmic for d = 3 and a — 0. 
As long as Eq > T, t can be considered as an extra 
dimension and the gauge invariance ([9]) prohibits gener- 
ating a relevant term oc <i>^. The other possible terms 
allowed by symmetry are irrelevant, i.e. the theory is re- 
normalizable. For e — I, Eqs. (8c I and (8e) still contain 



all the relevant terms. Moreover, term oc gi/c describes 
the long range Coulomb interaction and cannot be re- 
normalized. The term oc (72 is leading irrelevant and it 
describes the logarithmic interaction and its renormaliza- 
tion due to the virtual boson-antiboson pairs. 

We subdivide the fields A and $ into slow A, $ and 
fast Aq, $0 parts (In the first loop approximation cut-off 
procedure can be rather arbitrary, we treat energy Eq as 



a running cut-off), integrate in Eq. (8a) over Aq, $0 mak- 



ing a cumulant expansion in A, $ up to the fourth order. 



The action is reproduced with the couplings running as 
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(<7i + 92) 91 + 92 



52 



-52 D 
(the running of / is of no consequence) . Here 

/ = In {B/Eo) , 



(10a) 
(10b) 
(10c) 
(lOd) 
(lOe) 

(11) 



and Eqs. (8c), i8a are the initial conditions for Eqs. (^10 



The terms in curly brackets correspond to the /3-function 
for XY model in li + 1 dimension. 
RG flow should be stopped at ^ > Z*, 

h^ma^ik-Ja); |a(/a)|^l, h^HB/T), (12) 

where the quantum fluctuations loose importance either 
due to the finite temperature {l^, = It) or to depar- 
ture from the phase transition line (l^ = /„). In the 
former case the quantum RG analysis (in d + I dimen- 
sions) should be supplemented by the analysis of the d- 
dimcnsional classical fluctuations. 

This analysis deserves some discussion. Let us obtain 



the free energy from Eqs. (8b) and (8el by including only 
time independent fields: 

F= f d'^r{4 |VA|' + arT^ | A|' + brTT^c^T' \Af } 

(13) 



+ / d\ lA 



<i„IA|2^2 



d'^k |fc|''-l|$(fc)|2 

(27rp Att^ct 



5l,T + 52, T 



if) 



where the subscript T means that the couplings are calcu- 
lated at / = Zt- To obtain the canonical form one has to 
integrate e"^/"^ over all the configuration of <I>. Such in- 
tegration would immediately produce the term ~ — |A|^ 
resulting in the first order phase transition. Similar effect 
of the fiuctuating magnetic field was studied long ago[n]. 
However, the electrostatic potential is very different 
from the vector potential of the magnetic field. The 
gauge invariance prohibits the screening of the static vec- 
tor potential by the fluctuations of A. On the contrary, 
the static electrostatic potential can be screened [gauge 
invariance ([9]) does not allow to remove zero Matsubara 
component of $(?').] Explicit calculation of the static 
polarization operator adds extra term 



F„ 



6T 



l$* 



d'^k 

(^ 



(14) 



to the free energy F [this term is an effect of finite r 
and does not appear in RG analysis (10)]. As the result, 



the long-range fluctuations of ^ are massive, the last two 



lines in Eq. ( 18 ) can be neglected and we are left with 



a classical XF-model in d dimensions. The phase tran- 
sition is of the second order in rf = 3 and Bcrezinskii- 
Kosterlitz-Thouless in d — 2. Transition temperature 
estimated as 



Tr = Be" 



ailc)^-bih); 



(15) 



is shown by solid line on Fig. [T] 

First order phase transition at d — S- Second order 



phase transition (151 implies 5 > 0, so that the energy 
Ut{\A\) ^ a{le)£'\A\^ + 7T'c{krb{k)\A\\ (16) 



with £ = max(r, |A| 



l£ — In Y has a single minimum. 



In what follows, we solve the RG equations ( 10 1 and 



show that at some value oi It, b{lt) = 0- It means, that 
at T < Tt = B exp{—lt) the additional stable minimum 



appear in Eq. (16); the transition is of the first order 



theory is massive and renormalization terminates. The 
vicinity of Tt can be analyzed in the spirit of Ref. |10j . 
For d = 3 only 5 = (?i -I- (72 is important and we find 



'(0 = |, 5(/) = i^(^ln|, 



e = ^ + ^, (17) 



where gf~'^,b are given by Eqs. ( 8c|8f ), and ^* is found 
from 6gf=%/5 = F[\n{6gf=y5C)]- Function 



T{x) 



22 
25 



25 



cot 



/239 

To"' 



(18) 



changes sign at a; ss —0.62 .... Therefore, b changes the 
sign at finite It independently on the initial conditions, 
see Fig. nk. As the initial couplings, gf^^ ^ l, b < I 
the region of the possible phase coexistence, see Fig. [T] 
occupies significant part of the phase diagrams. 

First order phase transition at d = 2- The extrapo- 
lation of RG equations to e = 1 shows that for a weak 
Coulomb interaction, gf^^ ^ 1, RG flow has a stable 
fixed point, i.e. the SIT is of the second order. Indeed, 
for this case gi{l) = 21/3 + const, g2 — >■ e/6. Then, b 
tends to its fixed point value as & « 4e/5 -I- 0{l/gi) and 
the effect of the Coulomb interaction vanishes at I —>■ 00. 
Therefore, the effect of the Coulomb interaction may lead 
only to the logarithmic corrections to the usual power 
laws of d = 3 classical XY model. 

Situation changes in the opposite limit when the mu- 
tual capacitance significantly exceeds the intergrain dis- 
tance i.e. 52 ^ 51 ■ We can set in Eqs. (10) 51 = 0, 



then 32 rapidly reaches its fixed point 5/6e and for e = 1 
evolution of b{l) is governed by function (lls]): 



(19) 
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FIG. 1. Phase diagram of the JJA. The sohd and double- 
lines are the second and first order respectively. Insets: (a) 
the two parameter RG flow for 3D JJA; (b) two parameter 
RG flow for 2D arrays with logarithmic interaction; (c) the 
projection of three parameter RG flow for the logarithmic and 
"moderate" Coulomb interaction gi{l = 0) = 0.65. 



where b.^, ^ bg2. Accordingly, b always changes sign and 
the SIT is of the first order at T < Tt ~ E^oe"''*', see 
Fig.[l]D. If the Coulomb interaction is moderate, 32 — <?i, 
both kind of phase transitions are possible, see Fig. [l]:. 

Physical interpretation - The stability of both insulat- 
ing and superconducting states at T < Tt appears due 
to the competition of the effects of the long range inter- 
action of excitations against their Bose statistics. When 
the former is strong enough there exists a state formed 
by boson/antiboson dipoles. This stable state competes 
with the formation of the uniform Bose-condensate. 

In conclusion we have derived an effective field theory 
describing the SIT transition in granular superconductors 
and JJA. The RG analysis of this model demonstrated 
that the SIT is inevitably of first order at low temper- 
atures for all 3-dimensional and realistic 2-dimensional 
JJA. This may have very important experimental con- 
sequences. In particular, one can see hysteresis when 
changing e.g. magnetic field. The insulating and su- 
perconducting state can coexist and phase separation in 
space, especially in the presence of disorder, is possible. 
All these interesting phenomena deserve a separate study. 
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